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ABSTRACT 



In the first two papers of this series, we determined the position of Sgr A* 
on infrared images, by aligning the positions of red giant stars with positions 
measured at radio wavelengths for their circumstellar SiO masers. In this paper, 
we report detections of 5 new stellar SiO masers within 50" (2 pc) of Sgr A* and 
new and/or improved positions and proper motions of 15 stellar SiO masers. The 
current accuracies are ~ 1 mas in position and ~ 0.3 mas y" 1 in proper motion. 
We find that the proper motion of the central stellar cluster with respect to 
Sgr A* is less than 45 km s _1 . One star, IRS 9, has a three-dimensional speed 
of ~ 370 km s -1 at a projected distance of 0.33 pc from Sgr A*. If IRS 9 is 
bound to the inner parsec, this requires an enclosed mass that exceeds current 
estimates of the sum of the mass of Sgr A* and luminous stars in the stellar cusp 
by « 0.8 x 10 6 M . Possible explanations include i) that IRS 9 is not bound to 
the central parsec and has "fallen" from a radius greater than 9 pc, ii) that a 
cluster of dark stellar remnants accounts for some of the excess mass, and/or iii) 
that Rq is considerably greater than 8 kpc. 
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Subject headings: Galaxy: center - astrometry - masers - stars: dark matter: 
AGB and post AGB, variables 



Introduction 



Sagittarius A* (Sgr A*), the compact radio source at the center of the Galaxy (IBalick fc Brown 



19741 ). is almost surely a super- massive black hole. Infrared-bright stars orbit about the radio 



posit i on of Sgr A* and r equire ~ 3.9 x 10 6 M within a radius of ~ 50 AU (jEisenhauer et al. 



20051 : iGhez et al.l 120051 ). While orbiting stars move at thousands of km s , Sgr A* is 



essentially stationar y at the dynamical center of the Galaxy, moving < 2 km s 1 out of 
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2005), contains a significant fraction (> 10%) of the gravitational mass 



The Galactic Center stellar cluster contains red giant stars that are both strong radio 
sources (from circumstellar SiO maser emission) and bright infrared sources. Because these 
stars are visible at both radio and infrared wavelengths, one can use their radio positions, 
measured with respect to Sgr A*, to determine accurate scale, rotation, and distortion 
corrections for an infrared image. This allows the highly accurate radio reference frame to 
be transfered to the infrared images, improvi ng the quality o f the infra red astrome t ry an d 
precisely locating the position of Sgr A*. In IMenten et al. I (119971 ) and iReid et al.l (120031 ). 
hereafter Papers I & II, we developed this technique and determined the position of Sgr A* 
on diffraction-limited 2.2 fim wavelength images of the Galactic Center with an accuracy 
of ~ 15 mas. Locating Sgr A* on infrared images has been important for determining its 
emission in the presence of confusing stellar sources and verifying that foci of stellar orbits 
coincide with the position of Sgr A*. This links the radiative (compact radio) source with 
the gravitational source. 

We present new VLA observations of stellar SiO masers in the central cluster, updating 
their positions and, by more than doubling the observing time-span, significantly improving 
their proper motion determinations. We place the proper motions of infrared stars in the 
central cluster in a reference frame tied to Sgr A*. In §2 we describe the radio measurements 
of the positions and proper motions of 15 SiO maser stars, and in §3 we present the latest 
infrared positions and proper motions of those stars within 20" of Sgr A*. We use these data 
to transfer the infrared proper motions to a reference frame tied to Sgr A* in §4. Finally, 
in §5, we use the three-dimensional speeds and projected offsets of stars from Sgr A* to 
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constrain the combined mass of Sgr A* and the central (luminous and dark) stellar cluster. 



2. Radio Observations 

Over the period 1995 to 2006 we have searched for and mapped SiO masers associated 
with late-type stars that are projected near Sgr A*. We used the NRAO0 VLBA and VLA 
to measure accurately the relative positions of SiO maser stars and Sgr A*. 

Red giant and super-giant stars of late-M spectral class often exhibit SiO masers in 
their extended atmospheres. These masers are strongly variable ove r time scales of ~ 1 y. 



SiO m aser emission emanates from a typical radius of ~ 4 AU (e.g., iDiamond fc Kemball 



( 120031 ) ), which cor responds to ~ 0.5 mas at the assumed 8.0 kpc distance of the Galactic 



Center (jReidlll993l ). As our measurements were made over a time-span much longer than the 
stellar cycle, we cannot track individual maser features, and we accept an intrinsic stellar 
position uncertainty of about ±0.5 mas owing to possible variations across the maser shell. 
For a late-type super-giant, the stellar radius is considerably larger than for a red giant of 
similar spectral class, and its SiO masers are found at radii roughly an order of magnitude 
larger. Variation of the SiO masers in a super-giant can considerably degrade inferred stellar 
position and proper motion measurements (see discussion of IRS 7 in §2.1). 



2.1. VLA Observations 

Our VLA observations were conducted in the A- configuration in 1995 June (reported 
in Paper I), 1998 May and 2000 October/November (reported in Paper II) and 2006 March 
(reported here). Near the Galactic Center, SiO masers are likely detectable over a very 
wide range of velocities, probably exceeding 700 km s _1 . However, wide-band observations 
at the VLA are currently severely limited by the correlator. In order to obtain adequate 
spectral resolution and sensitivity, we chose to limit our velocity coverage and observe in 
seven 6.25 MHz bands (each covering rs 40 km s _1 excluding band edges). We observed 
in both right and left circular polarization for each band and obtained 64 spectral channels 
per band, resulting in uniformly weighted spectral channel spacings of about 98 kHz or 
0.67 km s _1 . 

Our 2006 observations were conducted on March 5, 18, & 19. The VLA had 24 antennas 



1 NRAO is a facility of the National Science Foundation operated under cooperative agreement by Asso- 
ciated Universities, Inc. 
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in operation, and the synthesized beam toward Sgr A* was about 86 x 33 mas elongated 
in the north-south direction. Except for occasional calibration sources, we pointed on the 
position of Sgr A*, allowing detection of masers within the primary beam of a VLA antenna 
(« 30" HWHM at 43 GHz). We observed by cycling among bands centered at LSR velocities 
of —346, —111, —73, —39, —1, +40 and + 75 km s _1 ; the latter six bands covered the LSR 
velocity range —131 to +95 km s _1 with only two small gaps. This setup allowed for deeper 
integrations for the previously known, or suspected, SiO maser stars than would be possible 
for a "wide-open" search. 

Initial calibration of the VLA data was done in a standard manner recommended by 
the Astronomical Image Processing System (AIPS) documentation. The flux density scale 
was based on observations of 3C 286, assuming 1.49 Jy for interferometer baselines shorter 
than 300 kA. Amplitude and bandpass calibration was accomplished with observations of 
NRAO 530, which had a flux density of 2.8 Jy during the 2006 observations. The visibility 
data were then self-calibrated (amplitude and phase) on Sgr A* for each individual 10- 
sec integration. This removed essentially all interferometer phase fluctuations, owing to 
propagation through the Earth's atmosphere, and placed the phase center at the position of 
Sgr A*. 

Sgr A* is known to vary occassionally in flux density by ~ 10% on hourly time scales 



(lYusef-Zadeh et al.ll2006l ). Were this to happen during our observations, self-calibrating on 
Sgr A* would introduce false gain (amplitude) variations of a similar magnitude and degrade 
somewhat the dynamic range of the SiO maser images. However, this is unlikely to shift the 
measured position of the masers significantly, as the positions are determined almost entirely 
by the phase data. 

We searched for new maser stars by making very large images, covering about ±50", or 
most of the primary beam of an individual VLA antenna at 43 GHz, about Sgr A*. This was 
done by limiting the range of uv-data used and resulted in maps with a nearly circular beam 
of 0.15". Typical rms noise levels in these images were near 9 mJy, allowing 6a detections 
of 54 mJy. Five new SiO maser stars were discovered: SiO-14, SiO-15, SiO-16, SiO-17, and 
IRS 19NW. 

Once the approximate location of a maser was known, either from previous observations 
or from the large images, we mapped each band with up to five small sub-images centered 
on known or suspected masers with emission in that band. These synthesized maps typically 
had single spectral-channel noise levels of about 5 mJy. We always included a sub-image 
for Sgr A* at the phase center of the interferometric data. By simultaneously imaging the 
stellar SiO masers and the continuum emission from Sgr A*, the strong continuum emission 
from Sgr A* did not degrade the detections of relatively weak SiO masers far from Sgr A*. A 
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composite spectrum from all seven observing bands for the 1998, 2000, and 2006 observations 
is shown in Fig. [TJ 

As described in Paper II, we obtained a single position for each star at each observing 
epoch by 1) fitting a 2-dimensional Gaussian brightness distribution to each spectral channel 
with detectable SiO emission, 2) averaging, using variance weighting, to obtain a best stellar 
position and estimated uncertainty (if the reduced x 2 was greater than unity, we increased 
the position uncertainties accordingly), and 3) correcting for differential aberration, an effect 
of < 1 mas for stars < 15 " of Sgr A*. 

We list the positions of the SiO maser stars, relative to Sgr A*, in Table [TJ We include 
the results from the 1995 VLA observations reported in Paper I, the 1998 and 2000 VLA 
observations reported in Paper II, and the 2006 results reported in this paper. Since stellar 
SiO masers are variable in strength over the period of the stellar pulsation and the sensitivity 
of each epoch's data differed somewhat, only the stronger sources are detected at all epochs. 

One of the stars, IRS 7, is a super-giant. As discussed in Paper II, it has SiO maser 
features spread over m 20 km s _1 and is highly variable. We would expect its stronger 
SiO maser peaks could be spread over a region of at least 10 mas. Thus, the positions 
determined from the brightest SiO maser feature (s) in IRS 7 may not indicate the stellar 
position to better than about 5 mas, and we have increased the position uncertainties for 
IRS 7 in Table [2] to allow for this possibility. 

We constructed spectra at the pixel of peak brightness for SiO masers detected in the 
1998, 2000, and 2006 VLA observations. These spectra are displayed in Figs. [2]&[3j Most 
of these SiO spectra are as expected for Mira variables located at the distance of the the 
Galactic Center: they show flux densities ^ 1 Jy covering a velocity range of 5 to 10 km s _1 
and strong variability over timescales of years. Fig. H] shows the Galactic Center region, with 
the positions and proper motions of the nine SiO maser stars that are projected within the 
inner 21" displayed. 

2.2. SiO Maser Proper Motions 

We determined stellar proper motions by fitting a variance-weighted straight line to the 
positions as a function of time from all of the available data compiled in Table [TJ These 
proper motion fits are given in Tableland displayed graphically in Figs.[2]&[3j The reference 
epoch for the proper motion solution was chosen as the variance-weighted mean epoch for 
each star, in order to obtain uncorrelated position and motion estimates. Since the estimated 
uncertainties for individual east-west and north-south positions were neither identical, nor 
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exactly linearly related, we chose a single, average reference epoch for each star (instead of 
a separate reference epoch for the east- west and north-south directions), which resulted in 
slight parameter correlations. 



Assuming a distance of 8.0 kpc to the Galactic Center (jReidlll993l ). we convert proper 
motions to linear velocities and, with the radial velocities, determine the full 3-dimensional 
speed, Vtotah °f eacn star with respect to Sgr A*. These speeds are given in Table [31 While 
the speeds on the plane of the sky are directly referenced to Sgr A*, the speeds along our 
line-of-sight are in the LSR reference frame. Thus, our values of V tota [ assume that Sgr A* 
has a near-zero line-of-sight speed with respect to the LSR. As no spectral lines have been 
detected from Sgr A*, there is no direct observational evidence supporting this assumption. 
However, the lack of a detectab le proper motion of Sgr A* suggests that it anchors the 
dynamical center of the Galaxy faeid Sz Brunthaler 2004 ) and should be nearly at rest in 
the LSR frame. 



Infrared Observations 



Near-IR observations were obtained on the 8.2-m UT4 (Yepun) of the ESO-VLT on 
Cerro Paranal, Chile, using the detector sy stem NAOS/CONICA (hereafter NACO) consist- 
ing of the adaptive optics system NAOS (IRousset et al. 1120031 ) and the 1024 x 1024-pixel 
near-IR camera CONICA Jflartung et al. Il2003h . We obtained 8 data sets in H and K bands 
with a scale of 27 mas/pixel (large scale) covering 5 epochs (May 2002, May 2003, June 2004, 
May 2005, April 2006). 

Each image covers a field of view (FOV) of 28 x 28 arcsec. During each observation 
the camera pointing was shifted so that the total FOV of a complete data set was between 
35 x 35 and 40 x 40 arcsec, centered on Sgr A*. In all cases the spatial resolution was (nearly) 
diffraction limited, leading to a typical FWHM of ~ 60 mas. Typical limiting magnitudes 
were 18 th mag at K band and 20 th mag at H band. All images were sky-subtracted, bad- 
pixel removed and flat-field corrected. In order to obtain the best signal-to-noise ratios and 
maximum FOV coverages in single maps, we combined all good-quality images belonging to 
the same data set to mosaics after correcti ng for instrumental geometric distortion. Details 
of the distortion correction will be given in lTrippe et al.l (lin preparation! ) . 



In order to establish an astrometric near-IR reference frame, we selected one high-quality 
mosaic (May 2005), extracted image positions for all detected stars and transformed them 
into astrometric coordinates using the positions of all 9 maser stars in the FOV interpolated 
to the epoch of the image. To compare positions among images of different epochs, we chose 
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an ensemble of ~ 600 "well-behaved" stars (i.e., stars that are well-separated from neighbors 
and are bright but not saturated) in the "master mosaic" and computed in each image all 
stellar positions relative to this ensemble. Due to varying FOVs, the number of stars usable 
for a given mosaic varied between about 400 to 600. 

In each image, positions were extracted by fitting stars with 2-dimensional elliptical 
Gaussian brightness distributions. Although over the entire FOV significant departures from 
isoplanicity occur, this effect elongates the stellar PSFs symmetrically and does not affect 
significantly the centroids of Gaussian-fitted positions. Proper motions were computed by 
fitting the positions versus times with straight lines. To assure statistically "clean" errors 
for the proper motions, outlier rejection and error rescaling were applied where possible. 
This led to typical position accuracies of ~ 1 mas and typical proper motion accuracies of 
« 0.3 mas y _1 . Unfortunately, this accuracy was not achieved for all maser stars, as some 
are very bright stars and are affected by detector non-linearity /saturation in some images; 
also the star most distant from Sgr A* (IRS19NW) was observed only in the last two epochs. 
Thus, some of the maser stars have errors larger than typical. 



4. Radio & Infrared Frame Alignments 

The proper motions of stars in the Galactic Center cluster, measured from infrared 
images, are relative motions only. One can add an arbitrary constant vector to all of the 
stellar proper motions without violating observational constraints. Until now, the "zero 
points" of the motions have been determined by assuming isotropy and removing the average 
motion of the entire sample. Since the radio proper motions are inherently in a reference 
frame tied directly to Sgr A*, one can use any one of the SiO stars, or the mean motion of 
a group of them, to place the infrared proper motions in Sgr A*'s frame. 

We have measured radio positions and proper motions, relative to Sgr A*, for the nine 
bright SiO maser stars that appear on the NACO images (ie, within « 20" of Sgr A*). The 
position and proper motion accuracies typically are ~ 1 mas and ~ 0.3 mas y _1 , respectively. 
This allows us to align the radio and infrared frames, both in position and in proper motion. 
Thus, not only can the location of Sgr A* can be accurately determined on infrared images, 
but also stellar proper motions from infrared data can be referenced directly to Sgr A*, 
without assumptions of isotropy or homogeneity of the stellar motions. 

The radio and infrared proper motions measured for the nine stars are listed in Table |H 
The nine stars have weighted mean differences (and standard errors of the means) of +0.66 ± 
0.21 mas y _1 toward the east, and —0.45 ±0.28 mas y _1 toward the north. These results are 
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qualitatively similar to those published in Paper II. Quantitatively, the differences between 
the IR motions of Paper II and this paper for some stars (notably IRS 12N) are greater than 
expected based on the quoted uncertainties. Since the IR motions in Paper II were based on 
2 epochs only, which did not allow for an internal check on the formal motion uncertainties, 
we believe those uncertainties were somewhat optimistic. 

Currently only one star, IRS 7, has a significant discrepancy between the radio and 
infrared motions in the both coordinates. This is the only super-giant star in the sample 
and, owing to its extreme brightness, the infrared measurements are compromised by detector 
saturation. Additionally, the radio measurements are subject to significant uncertainty from 
the large SiO maser shell size. After removing IRS 7, the weighted mean differences between 
the radio and IR motions change only slightly and become +0.63 ±0.21 mas y _1 toward the 
east, and —0.32 ± 0.18 mas y _1 toward the north. 

When comparing how well the IR frame matches the radio frame, we need to con- 
sider the statistical uncertainty of the average IR motion, which has been removed. For 
most epochs, the average motion is based on m 400 stars, each of which has a typical mo- 
tion of ~ 100 km s _1 . Thus, the mean IR motion should have an uncertainty of roughly 
100 km s -1 /\/400 ~ 5 km s -1 . Adopting the result with IR S 7 remove d, converting to 



linear speeds for a distance of 8.0 kpc to the Galactic Center (lReidlll993l ). and adding in 
quadrature a ~ 5 km s _1 uncertainty for the mean IR motion removed from each coordi- 
nate, implies that the infrared stellar cluster moves +24 ± 9 km s _1 toward the east, and 
— 12 + 9 km s _1 toward the north, with respect to Sgr A*. The northward component motion 
does not deviate significantly from zero; the eastward component formally presents a 2.7cr 
significance. Combining these components in quadrature formally yields a speed difference 
of 27 ± 9 km s _1 . However, at this time, we do not consider that we have firmly detected 
motion of the stellar cluster, and we adopt a 2a upper limit of 45 km s _1 for the proper 
motion of the stellar cusp with respect to Sgr A*. 



5. Enclosed Mass versus Radius from Sgr A* 

Estimates of the enclosed mass versus projected radius from Sgr A*, based on infrared 
stellar motions, rely on relative motions not tied directly to Sgr A*. Since, the 3-dimensional 
motions of the SiO masers in this paper are both very accurate and directly tied to Sgr A*, 
they provide valuable information on the enclosed mass within projected radii of 0.2 to 2 pc 
of Sgr A*. 

In Paper II, we derived a lower limit to the enclosed mass at the projected radius of each 
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star, assuming the stellar motions reflect gravitational orbits dominated by a central point 
mass. Given the 3-dimensional speed, V to t a h an d projected distance from Sgr A*, r pro j, for 
each star, we obtained a strict lower limit to the mass enclosed, M ench within the true radius, 
r, of that star from Sgr A*. For a given enclosed mass, semi-major axis and eccentricity 
(e), the greatest orbital speed occurs at pericenter for e pa 1. Since the projected pericenter 
distance cannot exceed the true distance, we obtained 



M 



end 



> 



V 2 r 

v _ total ' pro] 

2G 



Note that this enclosed mass limit is a factor of two lower than would be obtained for a 
circular orbit. This lower limit approaches an equality only when three criteria are met: 1) 
r pro j pa r, 2) the star is near pericenter, and 3) it has an eccentricity near unity. The a priori 
chance of any star satisfying all three of these criteria is small, especially since a star in a 
highly eccentric orbit spends very little time near pericenter. Thus, Eq. (1) provides a very 
conservative limit on enclosed mass. 

We evaluate the lower limit to M enc \ using Eq. 1 by adopting conservatively the smallest 
total velocity allowed by measurement uncertainties, i.e., by subtracting 2a from V to tai in Ta- 
ble |3] before calculating a mass limit. The mass limits, given in Table EH are mostly consistent 

with t he en closed mass versus projected dis t ance from Sgr A* given by lGenzel. Eckart. Ott. & Eisenhauer 
( 119971 ) and lGhez. Klein. Morris fc Becklinl (119981 ). For many of the stars, the lower mass lim- 
its are well below the estimated enclosed mass curves, as expected given the very conservative 
nature of the calculated limits. 



Our most significant lower mass limit is from IRS 9. In Paper II, we arrived at a 
limit > 4.5 x 10 6 M M , which exceeded the then favored model of a 2.6 x 10 6 black 
hole ( jGenzel. Eckart. Ott. fc Eisenhauer! 119971 : iGhez. Klein. Morris fe Becklinl Il998l ) , plus a 
0.4 x 10 6 M contribution from the central stellar cluster ( Genzel et al. 20031 ). by about 50%. 



With our improved proper motions, we now find a more stringent limit of > 5.1 x 10 6 M 
at a projected radius of 0.33 pc from Sgr A*. 



Fig. [5] displays our enclosed mass versus radius constraint based on the 3-dimension 
motion of IRS 9, along with other constraints in the recent literature. The current best 
estimate for the mass of the SMBH (Sgr A*) is (3.9 ± 0.2) x 10 6 M , for the distance 
to the Galactic Center, i?p, of 8.0 kpc. This m ass estimate comes fro m an unweighted 
average of the results of Eisenhauer et al.l (120051 ) and IGhez et al.l (120051 ) , based on stellar 
orbit determinations. Adding in a 0.4 x 10 6 M contribution fro m the central stellar cluster, 
based on a density profile of 1.2 x 10 6 (r/0.39 pc)" 14 M pc" 3 bv lGenzel et al. I (bo03h . yields 
4.3 x 10 6 M , still leaving a discrepancy of 0.8 x 10 6 M , for Rq = 8 kpc. Formally, this is 
about a 3a discrepancy, assuming an uncertainty of ±0.2 x 10 6 M in the mass estimate of 
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Sgr A* and an estimated ±30% uncertainty in the mass of the stellar cusp. 

Since we do not know the line-of-sight distance of IRS 9 from Sgr A*, one might be 
tempted to argue that r « ^T pro j and the star is simply sensing an enclosed mass of 1.0 x 
10 6 M from the central stellar cluster at that radius. However, the mass limit derived from 
Eq. (1) scales as r and would be approximately 10 7 M Q for r = 2r pro j = 0.66 pc. Thus, the 
mass discrepancy only increases for r > r pro j, as shown by the slanted line in Fig. [5] (but see 
§5.2). 

How can the lower limit to the enclosed mass provided by IRS 9 be explained? We now 
discuss some possibilities. 



5.1. Dark Matter in the Central Stellar Cluster 



One could explain the motion of IRS 9, were the central stellar cluster to cont ain dark 
matte r (in addition to Sgr A*) whose mass exceeds 0.8 x 10 6 M within r = 0.33 pc. iMorris 
(119931 ) estimates that ~ 10 6 M of "dark" stellar remnants (eg, white dwarfs, neutron stars 



black holes) cou l d hay e accumulated in the inner few tenths of a parsec of the Galaxy. 
Mouawad et al. I (120051 ) show that, with data available at the time, the orbital fit of star S2 



allows for (but does not require) 0.2 x 10 6 M of dark matter distributed within 0.001 pc of 
Sgr A*. Should such a dark component exist and extend to greater radii, it might explain 
some of the IRS 9 mass discrepancy. However, other estimates of th e total mass in black holes 
in th e central few tenths of a pc do not exceed ~ 0.2 x 10 6 M (IMiralda-Escude fc Gould 
2000l ; iFreitag. Amaro-Seoane fc Kalogera 1 120061 ) . Given these estimates and the evidence 
from other enclos ed mass indicators that do not support ~ 10 6 M of dark matter within 
w 1 pc of Sgr A* (IGenzel et al. 1120001 ). it seems unlikely that a dark component could explain 
more than a modest fraction of the IRS 9 mass discrepancy. 



5.2. IRS 9 not bound to the central parsec 

A critical assumption for calculating the minimum enclosed mass using IRS 9's space 
velocity (Eq. 1) is that it is in a bound orbit dominated by a central point mass. If IRS 9 is 
in a highly eccentric orbit with a semimajor axis greater than a few parsecs, this assumption 
can break down. In such a case, the star's space velocity could exceed the "local" escape 
velocity, based on the mass enclosed at its current radius, but still be bound at a larger 
radius. For example, a star could be bound by mass within « 10 pc of Sgr A*, but observed 
plunging into the inner few tenths of a parsec at a speed that makes it appear unbound. 
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This could happen via gravitational scattering starting either at sm all radii and increasing 



orbita l energy or at large radii and removing angular momentum. Ivan Langevelde et al. 



( 119921 ) suggest a similar explanation for three OH/IR stars with large radial velocities; these 
stars are seen projected tens of parsecs from Sgr A* and would require semimajor axes of a 
few kpc. 

Consider a Galactic Center star with a large semimajor axis and little angular momen- 
tum (e ~ 1), so that it essentially "falls" toward the Center. We derive the infall velocity 
for a star that starts "falling" at a radius, r max and reaches a radius of r . Assume a central 
point mass, Mbh, plus an extended component with density, p{r). The kinetic energy gained 
by a star falling from r to ro is equal to the difference in gravitational potential energy at 
those radii. Gravitational potential energy per unit mass, U m , for a spherically symmetric 
mass distribution has the properties that for mass interior to r 



and for mass exterior to r 



G f r 

U m = / p{r) Awr 2 dr 

r Jo 



tax r)( iy \ 

U m = -G I ^ 4vrr 2 dr 



r 



(For spherically symmetric systems, interior mass acts identically as a point mass equal to 
the enclosed mass located at the center of the distribution, while exterior mass results in zero 
gravitational force and a constant gravitational potential dependent on its radial position, 
but independent of the position of a "test" mass.) 



Following iGenzel et al. I (120031 ). p(r) = po(r/r ) a , where p = 1-2 x 10 6 M© pc~ 3 , 
ro = 0.39 pc, and a = —2.0 for r > ro pc and a ~ —1.4 for r < ro pc. This leads to 
an enclosed stellar mass within r of Mo/1.6, where M = 47rp ^o = 0-9 x 10 6 M Q , plus a 
contribution of Mbh from Sgr A*. Adding the contributions to the potential from different 
mass components for a star at radius r (for r > ro) gives 

U m = M BH - - / p(r) Anr 2 dr - - p(r) 4vrr 2 dr - G ^ 4vrr 2 dr . (2) 

r r J r J rQ J r r 

The first three terms on the right hand side of Eq. (2) sum the effects of the mass components 
interior to r, and the fourth term is the contribution from the mass exterior to r. Evaluating 
Eq. (2) we find 



U m = --M BH - - -M (- - 1) - -M fln(w/r Q ) - m(r/r )) . (3) 
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For a star "falling" from r max to ro, the kinetic energy per unit mass gained is equal to the 
difference in potential energy per unit mass. From Eq. (3), we find 

l v 2 = U m (r max ) - U m (r ) . (4) 



2 

where 

U m {r max ) = -— (m bh + M (^ - 0.38)) , 

and 

U m (r ) = (m bh + M (0.62 + \n(r max /r )) 

Evaluating Eq. (4) for r = 0.39 pc, which is a reasonable value for the 3-dimensional 
radius of IRS 9, gives v > 370 km s -1 for an initial radius r max > 9 pc. Thus, if IRS 9 is in 
a highly eccentric orbit that takes it out to a radius of > 9 pc, it could achieve its very high 
observed 3-D velocity without violating the enclosed masses estimated by other methods. 

A priori it might seem very unlikely that even one of 15 stars with detectable SiO 
masers would have such an orbit and be observed near closest approach to Sgr A* (where it 
spends little time). However, it is beyond the scope of this paper to evaluate the likelihood, 
especially with the limited statistics available at this time. 



5.3. The distance to the Galactic Center (Ro) exceeds 9 kpc 

Were the mass of Sgr A* > 4.7 x 10 6 M Q , no mass discrepancy would exist. The best 
current mass estimates are based on fitting orbits for many stars and should be robust. 
However, the greatest uncertainty in the mass of Sgr A* comes its strong d ependence on the 



adopt ed value of Rq = 8.0 kpc for the distance to the Galactic Center. lEisenhauer et al. 



(120051 ) derive central masses from orbit fitting of 4.06 x 10 6 M Q when adopting R = 8.0 kpc 
and 3.61 x 10 6 M for a best fit Ro = 7.62 kpc. These values suggest an enclosed mass 
M enc i oc -Ro' 4 - O ur mass limit based on IRS 9's 3-D motion would also increase with Ro, 
but more weakly. Since the LSR velocity is the dominant component in the 3-dimensional 
motion for IRS 9 (and is not dependent on Rq), our minimum mass estimate (Eq. 1) scales 
approximately as M end oc R^ 3 , mostly through r pro j. Allowing _R to increase to about 9 kpc 



removes the mass discrepanc y. However, such a large value for Ro seems very unlikely (IReid 
19931 : lEisenhauer et al.ll2005l ). 
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5.4. Non-zero VLsr f° r Sgr A* 



Were Sgr A* moving toward the Sun along the line-of-sight with a speed > 30 km s^ 1 , 
this would lower V to tai and, hence, the M enc i limit to > 4.4 x 10 6 . W hile it seems ver y 
unlikely that a super-massive object would have such a large motion (jReid et al. 2003 ). 
we now consider this possibility. One method to approach this problem is to average the 
velocities of large samples of stars very close to Sgr A*, assuming that this average would 
apply to Sgr A*. 

Our sample of SiO maser stars, which should be nearly complete in the LSR velocity 
range —131 to +95 km s _1 , does not show any obvious bias. Infrared observations of CO 
band-head velocities from late-type stars in the central parsecs yield average LSR velocities 
that are not statistically different from zero. For example, t he int egrated CO-band head 
velocities (within a 20" diameter aperture) of iMcGinn et al. I (119891 ) indicate positive (neg- 
ative) velocities at positive (negative) Galactic longitude, consistent with the direction of 
Galactic rotation, and a value of —10 ± 25 km s~ x at the position of Sgr A*. (However, 
these authors find possibly significant stellar velocities of —47 ± 8 km s -1 for four point- 
ing offsets perp e ndicu lar to the Galactic plane.) Individual stellar velocities compiled by 
Rieke fc Rieke I (119881 ) of 54 stars pr ojected within 6 pc o f Sgr A* have a mean velocity 
of —20 ±11 km s -1 . Alternatively, IWinnberg et al. I (119851 ) and ISjouwerman et al.l (119981 ) 
measured velocities of OH masers for 33 and 229 OH/IR stars, respectively, within about 
~ 40 pc of Sgr A*, which yield average velocities of +7 ± 11 and +4 ± 5 km s^ 1 . Overall, 
it appears that radial velocities of stars near Sgr A* suggest an average LSR velocity near 
zero, within « 20 km s _1 . 



5.5. IRS 9 is (or was) in a binary 

Were IRS 9 in a tight, massive binary, perhaps a significant portion of its space velocity 
might be contributed by internal orbital motion, possibly reducing its speed with respect 
to Sgr A*. However, we have observed IRS 9 for about 8 years and see no changes in its 
radial or proper motion velocity components. The spectra of IRS 9 shown in Fig. [2] are 
characteristic of Mira variables, which show variable emission over a range of 5 to 10 km s~ x 
about the stellar velocity. We estimate that the stellar radial velocity of IRS 9 has changed 
by less than 2 km s _1 over 8 years. Also, the proper motion components are well- fit by 
constant velocities, with 2a upper limits to accelerations of 0.4 and 0.6 mas y -2 (15 and 
23 km s _1 y _1 at 8.0 kpc) for the eastward and northward components, respectively. 

The observed changes in radial velocity (or proper motion components) would have 
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different characteristics depending on the relative values of the time span of the observations, 
At, and the orbital period, P. For At > P/4, we would have sampled large changes in the 
orbital mean anomaly and hence would have seen quasi-random changes of magnitude equal 
to the orbital speed. Our limit of < 2 km s" 1 change in the radial velocity over 8 y, would 
place a limit of ~ 2 km s -1 for the radial component of any orbital velocity. Alternatively, 
if At < P/4, then we could be sampling only a small portion of an orbit and detecting a 
velocity change might be difficult. However, setting At = 8 y requires P > 32 y. Since 
any bright companion for IRS 9 would have been observed, we adopt companion mass m 
of < 10 M . This would even allow for most black hole companions. For m < 10 M and 
P > 32 y, we find an upper limit for an orbital speed of < 20 km s _1 . Thus, a binary 
orbital contribution to the observed space velocity of IRS 9 could not exceed ~ 20 km s _1 
and likely would be considerably less. Thus, it is highly unlikely that the extreme velocity 
of IRS 9 could be explained as owing to a binary orbit. 

Could IRS 9 have been in a binary system and come unbound (or bound at a much larger 
radius as discussed in §5.2) after a close encounter with Sgr A*? A small number of "hyper- 
velocity" stars are thought to have been ejected from the Galactic Center in this manner 



(jHills Ill988l ; lYu fc Tremaine 1120031 ; iBrown et al. II2005I ). However, these are estimated to be 



very rare events (< 1 in 10 5 y) and we are statistically unlikely to be witnessing a newly 
created hyper- velocity star so close to Sgr A*. All hyper- velocity stars discovered to date are 
early-type main-sequence stars; they are found in the outer Galaxy and are moving at speeds 
of ~ 1000 km s _1 , even after climbing out of the gravitational potential of the inner Galaxy. 
Main sequence stars can survive the strong tidal forces experienced during close encounters 
with Sgr A*. However, IRS 9 is an AGB star and, thus, is a very extended (~ 300 solar 
radii) and low surface-gravity object. It is unclear if such a star could survive the ejection 
event, without losing its extended atmosphere. 



6. Conclusions 

We have now measured the radio positions and proper motions of 15 late type stars with 
SiO maser emission in the Galactic Center stellar cluster. All but two of these stars have been 
detected at three or more epochs and have measurement accuracies of ~ 1 mas in position 
and ~ 0.3 mas y _1 in proper motion. Nine of these stars have multi-epoch measurements 
of proper motions at infrared wavelengths. A comparison of the radio motions, which are 
relative to Sgr A*, with the infrared motions indicates that the stellar cusp moves with 
Sgr A* to within 45 km s -1 . 

The three-dimensional speeds and projected distances of individual stars from Sgr A* 
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yield lower limits to the enclosed mass. The enclosed mass limit for one star, IRS 9, exceeds 
current estimates of the combined mass of Sgr A* and the luminous stars in the cusp within 
the central parsec. This result is puzzling, but might be explained, for example, by a 
combination of i) a population of dark stellar remnants in the central parsec, ii) IRS 9 being 
on a plunging "orbit" with a semimajor axis 3> 1 pc, and/or iii) a value of Rq > 8 kpc. 



We thank S. Gillessen for comments on the paper, and A. Loeb and R. O'Leary for 
discussions on the extreme motion of IRS 9. 
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Fig. 1. — Composite spectrum of stellar SiO masers detected with the VLA in 1998 (blue), 
2000 (green) and 2006 (red) observations. Stars are identified at the top of the spectrum at 
their approximate stellar velocities. 
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Fig. 2. — Eastward (left panels) and northward (middle panels) position offsets from Sgr A* 
versus time for the eight SiO maser stars with negative LSR velocities. Solid lines are 
variance-weighted best-fit proper motions. The linear speed is indicated in each frame, 
assuming a distance of 8.0 kpc. Spectra from 1998 (blue), 2000 (green) and 2006 (red) 
observations are also shown (right panels). Star names are indicated in the upper right 
corner of the right panels. 
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Fig. 3. — Eastward (left panels) and northward (middle panels) position offsets from Sgr A* 
versus time for the seven SiO maser stars with positive LSR velocities. Solid lines are 
variance-weighted best-fit proper motions. The linear speed is indicated in each frame, 
assuming a distance of 8.0 kpc. Spectra from 1998 (blue), 2000 (green) and 2006 (red) 
observations are also shown (right panels). Star names are indicated in the upper right 
corner of the right panels. 
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Fig. 4. — Infrared (K-band) image of the central ±20" of the Galactic Center taken in 2005, 
with east to the left and north to the top. SiO maser stars within this region are circled 
and their proper motions relative to Sgr A* are indicated with arrows. The vertical bar and 
arrow at the right of the image indicate the linear and motion scales for Rq = 8.0 kpc. The 
location of Sgr A* is indicated by the circle at the center of the image. 
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Fig. 5. — Constraints on the enclosed mass as a function of r adius (R) from Sgr A*. The 
estimates labeled "E" and "G" are from fitting stellar orbits by lEisenhauer et al.l (120051 ) and 
Ghez et al.l (120051 ) . respectively, a nd the estimate labe l ed "B " is from a statistical analysis 
of the "clockwise stellar disk" by iBeloborodov et al. I (120061 ). Our lower limit on enclosed 
mass from the 3-dimensional motion of IRS 9, assuming that it is bound in the region 
dominated gravitationally by Sgr A* (see §5.2), is indicated by a sloping line and arrow. 
The minimum distance of IRS 9 from Sgr A* is its projected distance of 0.33 pc. The 
sloping line corresponds to the mass limit for reasonable values of the unknown line-of-sight 
distance of IRS 9 from Sgr A*; this increases both the radius and the mass limit linearly. 
The dashed line ind icates the combined contribution of a point mass of 3.9 x 10 6 M© and a 
central stellar cusp (iGenzel et al. II2003I ). The uncertainty in the dashed line is dominated by 
the uncertainty of ±0.2 x 10 6 M in the point-mass, with a smaller contribution of perhaps 
±0.1 x 10 6 M Q from the central stellar cusp at R = 0.33 pc. All values assume Rq = 8.0 kpc; 
the discrepancy between the dashed line and the IRS 9 limit cannot be removed for values 
of R < 9 kpc. 
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Table 1. SiO Maser Astrometry 



Star 


Vlsr 




ABy 


Epoch 


Telescope 




( km s- 1 ) 


(arcsec) 


(arcsec) 


(y) 




IRS 9 


-342 


5.6515 ± 0.0048 


-6.3589 ± 0.0060 


1998.39 


VLA 






5.6501 ± 0.0007 


-6.3509 ±0.0013 


1998.41 


VLA 






5.6589 ±0.0011 


-6.3454 ±0.0017 


2000.85 


VLA 






5.6742 ± 0.0004 


-6.3347 ± 0.0009 


2006.20 


VLA 


IRS 7 


-114 


0.0403 ± 0.0080 


5.5829 ±0.0130 


1995.49 


VLA 






0.0387 ± 0.0023 


5.5676 ± 0.0049 


1998.39 


VLA 






0.0378 ± 0.0043 


5.5495 ±0.0014 


1998.41 


VLA 






0.0342 ± 0.0016 


5.5414 ± 0.0030 


2000.85 


VLA 






0.0326 ± 0.0007 


5.5238 ± 0.0013 


2006.20 


VLA 


SiO-14 


-112 


-7.6648 ± 0.0012 


-28.4528 ± 0.0020 


1998.41 


VLA 






-7.6596 ± 0.0005 


-28.4530 ± 0.0008 


2000.85 


VLA 






-7.6485 ± 0.0005 


-28.4586 ± 0.0009 


2006.20 


VLA 


IRS 12N .... 


- 63 


-3.2519 ± 0.0005 


-6.8814 ± 0.0005 


1996.41 


VLBA 






-3.2541 ± 0.0005 


-6.8877 ±0.0006 


1998.39 


VLA 






-3.2543 ± 0.0010 


-6.8860 ± 0.0011 


1998.41 


VLA 






-3.2554 ± 0.0009 


-6.8936 ± 0.0012 


2000.85 


VLA 






-3.2626 ± 0.0009 


-6.9073 ± 0.0019 


2006.20 


VLA 


IRS 28 


- 55 


10.4702 ± 0.0030 


-5.7884 ± 0.0050 


1998.41 


VLA 






10.4693 ± 0.0010 


-5.7956 ± 0.0024 


2000.85 


VLA 






10.4809 ± 0.0007 


-5.8254 ± 0.0010 


2006.20 


VLA 


SiO-15 


- 36 


-12.4253 ±0.0023 


-11.0794 ±0.0054 


2000.85 


VLA 






-12.4385 ±0.0012 


-11.0753 ±0.0015 


2006.20 


VLA 


IRS 10EE ... 


- 27 


7.6821 ± 0.0030 


4.2125 ±0.0050 


1995.49 


VLA 






7.6837 ± 0.0005 


4.2194 ± 0.0005 


1996.41 


VLBA 






7.6841 ± 0.0005 


4.2146 ± 0.0009 


1998.39 


VLA 






7.6837 ± 0.0005 


4.2157 ±0.0005 


1998.41 


VLA 






7.6845 ± 0.0005 


4.2099 ± 0.0005 


2000.85 


VLA 






7.6840 ± 0.0005 


4.1990 ±0.0005 


2006.20 


VLA 


IRS 15NE ... 


- 12 


1.2256 ±0.0140 


11.3108 ±0.0210 


1995.49 


VLA 






1.2302 ±0.0005 


11.3315 ±0.0005 


1996.41 


VLBA 






1.2249 ±0.0017 


11.3193 ±0.0019 


1998.39 


VLA 






1.2270 ±0.0005 


11.3201 ±0.0006 


1998.41 


VLA 






1.2228 ±0.0011 


11.3024 ±0.0025 


2000.85 


VLA 






1.2112 ±0.0005 


11.2761 ±0.0010 


2006.20 


VLA 


SiO-16 


+ 8 


-26.4237 ± 0.0020 


-34.4520 ± 0.0027 


1998.41 


VLA 






-26.4216 ±0.0006 


-34.4412 ±0.0011 


2000.85 


VLA 
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Table 1 — Continued 



Star 


Vlsr 
( km s- 1 ) 


AO x 
(arcsec) 


ABy 

(arcsec) 


Epoch 

(y) 


Telescope 






if? Him _i_ n nnnc 
— 26.4191 ± U.UUUb 


—34.4553 ± 0.0013 


2006.20 


\ n A 
V i^A 


SiO-6 


+ 52 


35.1982 ±0.0090 
35.2207 ± 0.0029 
35.2323 ± 0.0006 
35.2451 ±0.0010 


30.6567 ± 0.0140 
30.6537 ± 0.0050 
30.6593 ± 0.0010 
30.6702 ± 0.0028 


1995.49 
1998.39 
2000.85 
2006.20 


VLA 
VLA 
VLA 
VLA 


SiO-17 


+ 53 


8.0427 ± 0.0005 
8.0624 ± 0.0005 


-27.7034 ±0.0014 
-27.6852 ±0.0009 


1998.41 
2006.20 


VLA 
VLA 


SiO-11 .... 


+ 70 


1.7121 ±0.0040 
1.7379 ±0.0023 
1.7401 ±0.0005 
1.7462 ±0.0005 


40.2614 ± 0.0060 
40.2681 ± 0.0032 
40.2794 ±0.0011 
40.2914 ± 0.0006 


1995.49 
1998.39 
2000.85 
2006.20 


VLA 
VLA 
VLA 
VLA 


IRS 17 


+ 73 


13.1501 ±0.0026 
13.1415 ±0.0013 


5.5651 ±0.0025 
5.5611 ±0.0021 


2000.85 
2006.20 


VLA 
VLA 


SiO-12 .... 


+ 82 


-18.8645 ±0.0190 
-18.8235 ±0.0028 
-18.8212 ±0.0017 


42.4905 ± 0.0290 
42.4686 ± 0.0032 
42.4459 ± 0.0033 


1995.49 
2000.85 
2006.20 


VLA 
VLA 
VLA 


IRS 19NW 


+ 84 


14.5518 ±0.0011 
14.5532 ±0.0015 
14.5607 ± 0.0005 


-18.4619 ±0.0012 
-18.4683 ±0.0031 
-18.4656 ±0.0009 


1998.41 
2000.85 
2006.20 


VLA 
VLA 
VLA 



Note. 



VLBA 



AQ X and AS y arc 
in the J2000 coordi 
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Table 2. SiO Maser Proper Motions 



Star 


Vlsr 


AO* 


ABy 






Epoch 


Number 




( km s- 1 ) 


(arcsec) 


(arcsec) 


(mas y _1 ) 


(mas y -1 ) 


(y) 


Epochs 


IRS 9 


-342 


+5.6655 ± 0.0003 


-6.3407 ± 0.0007 


+3.06 + 0.10 


+2.11 + 0.19 


2003.34 


3 


IRS 7 


-114 


+0.0336 ± 0.0050 


+5.5401 ± 0.0050 


-0.58 + 0.50 


-3.52 + 0.54 


2004.37 


4 


SiO-14 .... 


-112 


-7.6554 ± 0.0003 


-28.4553 ± 0.0006 


+2.08 + 0.12 


-0.94 + 0.20 


2002.89 


3 


IRS 12N .... 


- 63 


-3.2537 ± 0.0003 


-6.8864 ± 0.0003 


-1.06 + 0.10 


-2.70 + 0.17 


1998.17 


4 


IRS 28 


- 55 


+10.4784 + 0.0011 


-5.8190 + 0.0010 


+2.00 + 0.38 


-5.29 + 0.42 


2005.00 


3 


SiO-15 


- 36 


-12.4372 + 0.0011 


-11.0757 + 0.0015 


-2.47 + 0.98 


+0.77 + 2.10 


2005.68 


2 


IRS 10EE ... 


-27 


+7.6840 ± 0.0003 


+4.2114 ± 0.0003 


+0.04 + 0.07 


-2.09 + 0.07 


2000.24 


5 


IRS 15NE ... 


- 12 


+1.2257 ± 0.0003 


+11.3171 + 0.0004 


-1.96 + 0.07 


-5.68 + 0.12 


1998.92 


5 


SiO-16 


+ 8 


-26.4207 ± 0.0004 


-34.4478 ± 0.0043 


+0.49 + 0.15 


-1.84 + 1.52 


2002.82 


3 


SiO-6 


+ 52 


+35.2333 + 0.0011 


+30.6605 ± 0.0009 


+2.58 + 0.43 


+1.99 + 0.52 


2001.43 


4 


SiO-17 


+ 53 


+8.0560 ± 0.0004 


-27.6911 + 0.0008 


+2.53 + 0.18 


+2.34 + 0.42 


2003.68 


2 


SiO-11 .... 


+ 70 


+1.7441 + 0.0014 


+40.2871 ± 0.0007 


+1.30 + 0.46 


+2.45 + 0.25 


2004.38 


4 


IRS 17 


+ 73 


+13.1442 + 0.0012 


+5.5624 ± 0.0016 


-1.61 ± 1.08 


-0.75 ± 1.22 


2004.49 


2 


SiO-12 .... 


+ 82 


-18.8227 + 0.0030 


+42.4559 ± 0.0023 


+0.77 + 1.14 


-4.24 + 0.84 


2003.83 


3 


IRS 19NW .. 


+ 84 


+14.5578 ± 0.0005 


-18.4647 + 0.0012 


+1.19 + 0.14 


-0.43 + 0.31 


2003.79 


3 



Note. — For soi 
have been correctei 
proper motions, rel; 
by +0.010" northw 
increased to 0.005" 
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Table 3. 3-Dimensional Stellar Motions & Enclosed Mass Limits 



Star 


VLsr 


V x 


Vy 


Vtotal 








( km s- 1 ) 


( km s- 1 ) 


( km s- 1 ) 


( km s- 1 ) 


(PC) 


(10 6 M ) 


IRS 9 ... 


-342 ± 3 


116 ± 4 


80 ± 7 


370 ± 3 


0.33 


> 5.1 


IRS 7 ... 


-114 ±3 


-22 ± 19 


-133 ±20 


177 ± 16 


0.21 


> 0.5 


SiO-14 ... 


-112 ±3 


79 ± 5 


-36 ± 8 


142 ± 4 


1.14 


> 2.4 


IRS 12N ... 


-63 ±3 


-40 ± 4 


-102 ± 6 


127 ± 6 


0.30 


> 0.5 


IRS 28 ... 


-55 ±3 


76 ±14 


-201 ± 16 


221 ± 15 


0.46 


> 2.0 


SiO-15 ... 


-36 ±3 


-94 ± 37 


29 ±80 


105 ± 40 


0.65 


> 0.0 


IRS 10EE ... 


-27 ±3 


-5± 3 


-82 ± 3 


87 ± 3 


0.34 


> 0.3 


IRS 15NE ... 


-12 ± 3 


-74 ± 3 


-215 ± 5 


228 ± 4 


0.44 


> 2.5 


SiO-16 ... 


8 ± 3 


19 ± 6 


-70 ± 58 


73 ±55 


1.68 


> 0.0 


SiO-6 ... 


52 ±3 


98 ± 16 


75 ± 20 


134 ± 16 


1.81 


> 2.2 


SiO-17 ... 


53 ± 3 


96 ± 3 


89 ± 8 


141 ± 6 


1.12 


> 2.2 


SiO-11 ... 


70 ±3 


49 ± 17 


93 ± 9 


126 ± 10 


1.56 


> 2.1 


IRS 17 ... 


73 ± 3 


-61 ±41 


-28 ± 46 


99 ±29 


0.55 


> 0.1 


SiO-12 ... 


82 ± 3 


29 ±43 


-160 ±32 


183 ± 29 


1.80 


> 3.3 


IRS 19NW ... 


84 ±3 


45 ± 5 


-16 ±12 


97 ± 4 


0.91 


> 0.8 



Note. — V x and 
total speed of the s 
assume a distance t 
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Table 4. Radio-Infrared Proper Motion Alignment 



c tal . ..Radio ..Radio ,.IR ,,IR .Dif Dif 

otar jx x fi y \i x fi y [i x ^ty 

(mas y _1 ) (mas y -1 ) (mas y -1 ) (mas y _1 ) (mas y _1 ) (mas y _1 ) 



IRS9 ... 


3.06 ±0.10 


2.11 ±0.19 


4.44 ±0.55 


1.10 ±0.58 


1.38 ±0.56 


-1.01 ±0.61 


IRS7 ... 


-0.58 ±0.50 


-3.52 ±0.54 


1.86 ± 1.18 


-7.33 ± 0.72 


2.44 ± 1.28 


-3.81 ±0.90 


IRS12N ... 


-1.06 ±0.10 


-2.70 ±0.17 


-0.77 ±0.66 


-3.39 ± 0.39 


0.29 ±0.67 


-0.69 ±0.43 


IRS28 ... 


2.00 ±0.38 


-5.29 ±0.42 


2.27 ±0.35 


-5.85 ± 0.31 


0.27 ±0.52 


-0.56 ±0.52 


SiO-15 ... 


-2.47 ±0.98 


0.77 ±2.10 


-2.18 ±0.51 


-0.56 ± 0.12 


0.29 ± 1.10 


-1.33 ± 2.10 


IRS10EE ... 


0.04 ± 0.08 


-2.09 ±0.07 


0.73 ±0.23 


-1.92 ± 0.27 


0.69 ±0.24 


0.17 ±0.28 


IRS15NE ... 


-1.96 ±0.07 


-5.68 ±0.12 


-2.40 ±0.48 


-6.29 ± 0.35 


-0.44 ± 0.49 


-0.61 ±0.37 


IRS17 ... 


-1.61 ± 1.08 


-0.75 ± 1.22 


0.20 ±0.59 


-1.67 ±0.62 


1.81 ± 1.23 


-0.92 ± 1.37 


IRS19NW ... 


1.19 ±0.14 


-0.43 ± 0.31 


-0.60 ± 3.08 


-0.54 ± 3.47 


-1.79 ±3.08 


-0.11 ±3.48 



Note. — fx x and 
(infrared minus rad 
motions arc relative 



